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FIXED-DOMAIN ASYMPTOTICS FOR A SUBCLASS OF 
MATERN-TYPE GAUSSIAN RANDOM FIELDS 

By Wei-Liem Loh 

National University of Singapore 

Stein [Statist. Sci. 4 (1989) 432-433] proposed the Matern-type 
Gaussian random fields as a very flexible class of models for computer 
experiments. This article considers a subclass of these models that are 
exactly once mean square differentiable. In particular, the likelihood 
function is determined in closed form, and under mild conditions 
the sieve maximum likelihood estimators for the parameters of the 
covariance function are shown to be weakly consistent with respect 
to fixed-domain asymptotics. 

1. Introduction. In the modeling of computer experiments, it has be- 
come rather common practice to approximate the deterministic response as 
a realization of a stochastic process. In this regard, Sacks, Welch, Mitchell 
and Wynn [10] and Sacks, Schiller and Welch [9] proposed modeling using 
a Gaussian random field X(x), x G [0, l] d , with a multiplicative covariance 
function, 

d 

Cov(X(x),X(y)) = a 2 J] exp(-0 t |xt - ytf) 

(1) 

Vx = (xi, . . . , x d )', y=(y 1 ,..., y d )' G [0, If, 

where 7 G (0, 2] , 0\ , . . . , 9 d and a 2 are strictly positive parameters. Ying [16, 
17] investigated the fixed-domain asymptotic properties of the maximum 
likelihood estimators of the covariance function when 7=1. In particular, 
he proved that the estimators are strongly consistent and asymptotically 
normal under mild conditions, van der Vaart [13] showed that when 7 = 
1 and al = 2, the maximum likelihood estimators are also asymptotically 
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efficient. Recently, Loh and Lam [6] showed that, when 7 = 2, sieve maximum 
likelihood estimators of 8i,...,6 d are strongly consistent using a regular 
sampling grid. 

Stein [11] observed that the Gaussian model given by (1) may have some 
undesirable properties. In particular, for 7G (0,2), the Gaussian random 
field with covariance function as in (1) will not be mean square differen- 
tiable. However, for the case 7 = 2, it is infinitely mean square differen- 
tiable. Not allowing for processes that are differentiable but not infinitely 
differentiable may be unnecessarily restrictive. Stein then suggested using 
a Gaussian random field model, X(x), x S [0,l] d , with the multiplicative 
Matern-type covariance function 

(2) Co V (X(x),X(y)) = g 2a ^% 2)e 2 a (0t\x t - yt\TK a (9 t \x t - y t \) 

Vx = (xi,...,x d )',y=(y 1 ,...,y d )' £ [0,l] d , 

where a, (j), 8±, . . . , d are positive constants and K a is the modified Bessel 
function of the second kind (see, e.g., [3], pages 222-223). The interesting 
parameter is a, where X will be m times mean square differentiable if and 
only if a> m. 

Due to the dependence of the models with a multiplicative covariance 
function on the choice of coordinate axes, Stein ([12], pages 48-55) later ad- 
vocated the use of Gaussian random field models with isotropic Matern-type 
covariance functions. This class of covariance functions was first proposed 
by Matern in 1960 as a reasonable class of models for isotropic random fields 
(see [7]). Stein ([12], Section 6.7) investigated the performance of maximum 
likelihood estimators for the parameters of a periodic version of the Matern 
model with the hope that the large sample results for this periodic model 
will be similar to those for nonperiodic Matern-type Gaussian random fields 
under fixed-domain asymptotics. 

This article assumes that X(x),x E [0, l] d , is a mean-zero Gaussian ran- 
dom field with covariance function given as in (2), where a = 3/2 and 
<f>, 9x, . . . , 9 d are unknown positive constants. Since K 3 / 2 ( s ) = (vrs _3 /2) 1 / 2 (l + 
s)e~ s (see [15], page 747), we observe that the covariance function of X(x) 
is 

Cov(X(x),X(y)) = ^ nU + 0t\xt ~ iftDe-* 1 *-* 1 

(3) f d t=1 

yx = (xi, . . . , x d )', y=( yi ,..., y d )' G [0, l) d . 

Remark. For reasons of mathematical tractability, we have assumed 
a multiplicative Matern-type covariance function and not an isotropic one. 
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The class of multiplicative Matern-type covariance functions has nonethe- 
less been used in, for example, [14]. A referee noted that the multiplicative 
Matern-type covariance function and an isotropic one have a major difference 
under fixed-domain asymptotics. In particular, Zhang [18] recently proved 
that consistent estimators do not exist for all parameters in the latter case 
for d = 1, 2 or 3. 

We are concerned with the estimation of (f),0i, . ■ . ,9d using observations 
that are taken from the above random field on a regular grid, that is, 

(4) I x( -,...,-) :l<i t <n,l<t<d 

I \ n n J 

where re is a strictly positive integer. 

For simplicity, we order the elements of the set (4) lexicographically as 
an n d x 1 column vector X n . Thus, the element X(ii/n, . . . ,id/n) precedes 
the element X(ji/n, . . . ,jd/n) in X n if and only if there exists a 1 < k < d 
such that it = jt whenever 1 < t < k and it < jk- Then the covariance matrix 
^<f>,ei,...,6 d ;n of X n is given by 

d ±d d 
^<j>,9i,...,e d ;n = 2 d#3 ... ^3 09 R 8t,n, 

where the symbol "0" denotes the Kronecker product (see [2], page 599), 
and for each 1 < t < d Rg un denotes the n x n matrix whose (i,j)ih ele- 
ment is (1 + 9 t \i - j\/n)exp[-0 t \i - j\/n]. Since X n ~ N n d(0, '^ < j,,e 1 ,...,e d ;n)> 
the likelihood function is 

L n {^9 1 ,...M = {^r nd/2 \^eu...fi^ 

and the log-likelihood satisfies 



21ogL n (0,0i,...,0 d ) 

-n d log(27r) -n d log 



( 7r 



d±d 



(5) - ■>' 1 ^-" d ~^¥er^i 

d 

log (g)Rg t>n 
t=l 



-1 



\j4 d X'n\ 



TT" 



n ■ 



The rest of this article is organized as follows. Theorem 1 in Section 2 
gives an exact closed- form formula for the determinant of the matrix R$ tn . 
Similarly, Theorems 2-4 in Section 3 establish an exact closed-form formula 
for the inverse of Re t) n- It should be noted that the closed-form formulas 
(especially the inverse) are not simple but fortunately are amenable to the- 
oretical analysis with the help of a mathematical software system such as 
Mathematica [15] that has symbolic computation capability. 

In Section 4 we simplify the exact formulas of the previous two sections 
via asymptotic approximations. These approximations are very sharp in that 
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the error is of the order 0(c n ) for some constant < c < 1. In particular, 
Theorem 5 gives an asymptotic approximation of the determinant \Re t) n\ 
and Theorems 6 and 7 give asymptotic approximations to the elements of 
the inverse Ra „■ 

Using the results of Section 4 and Mathematica, Theorem 8 in Section 
5 computes a large sample approximation of the Fisher information matrix 
for the parameters (f>,9±, . . . ,9d- 

Section 6 gives a simple consistent maximum likelihood-type estimate 
e d ( see Theorem 9) for cj) even if the parameters 9i, . . . ,9^ are misspeci- 
fied. Section 7 defines a sieve maximum likehood estimator ((f), 9i, . . . , 9d) for 
(4>,9i, . . . , 9d)- Theorem 10 establishes the weak consistency of (0, 9\, ... , 9^) 
under mild conditions when d > 3. 

Appendix A contains technical results that are needed in the main body of 
this article. Many of the results found in Appendix A, though conceptually 
simple, involve extremely complicated computations and series expansions 
and the use of Mathematica is critical here. Appendix B contains the defi- 
nitions of fourteen constants that are used in the exact expression of Rq^ti 
in Section 3. 



2. Determinant. In this section our main aim is to evaluate the deter- 
minant of Rq n f° r some strictly positive constant 9. For simplicity, we write 
w = 9/n,u = e~ w , and for 1 < m < n Re t n:m denotes the m x m matrix 
whose (i, j')th element is {Re,n\m)i,j = (1 + \i — j\w)v} l ~^ , for all 1 < i,j < m. 
Clearly, we have Re, n] n = Re,n- 

Proposition 1. The determinant of Rg in - m is given by 

\Re,n-,i\ = 1, 

and for 2 < m < n it satisfies the recurrence relation 
\Re,n;m\ = (-l) m ~ 2 {[l + (m - 2)w] - (1 + w)[l + (m - l)w]u 2 }u m ~ 2 V™" 2 

m— 3 

,n;m— k— 1 \t~ 1 > 

fc=0 

where 

ro, ifk = -2, 

T k = I (w — l)u + (1 + w)u 3 , if k = -I, 

{ (1 + kw)u k - 2[1 + (k + l)w]u k+2 + [1 + (k + 2)w]u k+i , if k > 0. 

Proof. Using linear algebra elementary row operations, we observe 
that, for each 1 < m < n, Re^m can be reduced to the m x m almost upper 
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triangular matrix Ag tn - m , where 



i,n;m 



(1 + w)u 





(l + w)u (l + 2w)u 2 (l + 3w)u 3 
(l + w)u (l + 2w)u 2 





V o 



i 

T-l 





To 
T-l 







Tl 
Tl) 






(l + (m - rjwy 1 - 1 \ 



(1 + (m - 2)w)u m ' 2 

Tm-i 
I'm — 4 



n 

TO 



J 



Now use elementary row operations again to reduce Ag tn;m to a diagonal 
matrix and then use the fact that the determinant of a diagonal matrix is 
equal to the product of its diagonal elements. This proves Proposition 1. □ 

The following theorem gives a closed form expression for \Re.n-m\- 

Theorem 1. For 2<m<n, 
\Re ,n;m\ 

(aai + u 2 b) \-{b + r_i<u)ai] m - 2 - (tia 2 + u 2 b) [-(b + r_iu)a 2 ] m - 2 



where 



a\ — (%2 



a = 1 - 2u 2 + u 4 - 2u 2 w + 2u 4 w, 
b = 2u 2 — u 4 + w — u 4 w — 1, 
a = 1 — u 2 — 2u 2 w — u 2 w 2 , 
b = u 2 + w + u 2 w — 1, 



a 2 



-(a — 2t_iu) — yj(a — 2r_iu) 2 — 4t_iu(6 + r_ii 
2(6 + T-m) 



-(a — 2t_ 1 u) + ^(a — 2r_in) 2 — 4t_i«(& + t_ii 
2(6 + t_iu) 



Proof. It is convenient to define 

(6) T * = ( (1 + kw)u k - (1 + w)[l + (k + l)w]u k+2 , iik>0, 

{ ' Tfe \0, otherwise. 

Then Tq = \Re, n ;2\- We observe from Proposition 1 that, for 2 < m < n, 

m— 1 

(7) \Re,n;m\ = Tm-2( _r -l) m ~ 2 + X! r m-fc-l(-1"-l) m ~ _1 |-Re,n;fc| • 

k=2 



G 
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Applying (7) recursively to the right-hand side of itself, we obtain 

\Re,n;m\ = ^m-2 (~ T -1 J™' 2 + Ttu-^Tq (-T-1 ) m ~ 3 
m— 1 

+ E r m-fci-l(-T-l) m ~ fcl ~ 1 |-R(9,n;fcil 
fci =3 

^* ( r \m—2 i _ _ \m-3 

= r m-2(- r -i) +r m _ 3 r (-r_i) 

m— 1 

+ £ r^^-iC-r-O" 1 -* 1 - 1 

fci=3 

fci-l 

+ E 7 ^-*2-l(- r -l)* 1 " fca " 1 |^,n i fc a l[ 
fc 2 =3 J 

m— 1 

= {- T -l) T m-2 + (- T -l) T m-fci-l T fc 1 -2 

fci =2 

m— 1 fei— 1 

+ E E T m-k 1 -lT~k 1 -k 2 -l{— T -l) m ~ k2 ~ 2 \Re,n;k 2 \- 
fc 1= 3 fc 2 =2 

Now continuing this argument repeatedly, we obtain 

m— 1 

\Rd,n;m\ = (- r -l) m_2 7"m-2 + (-T-l)™" 3 ^ T m _ kl _iT^_ 2 

fei =2 

m— 1 fci— 1 

+ (-r„i) m ~ 4 5] 5] r m _ fcl _ir fcl _ fc2 _iTfc 2 _2 

fc 1= 3/c 2 =2 
m— 1 fci— 1 fc 2 — 1 

+ (- T -l) m_5 E E E Tm - fc l- irfc l- fc 2-l r fc2-fc3-lT'fc3-2 
fc 1= 4 fe 2 =3 fc 3 =2 

m— 1 fci— 1 — 1 

H h(-r_i) ^ zJ ' E r m-fci-l7fci-fc 2 -l--' 

fci=m— 2 fc2=™- — 3 fc m _3=2 

X T fcm-4-fcm-3-l T fc m _3-2 

+ r m - 2 r *. 

(8) 

Hence, it follows from (8) and Lemma A.l (see Appendix A) that \Re, n ;m\ 
equals the coefficient of z m ~ 2 in the series expansion of 

(-T„ lU ) m - 2 G T *(z) + (-T^u) m ^zG T (z)G T *{z) + (-T„ lU ) m - A z 2 G 2 T (z)G T *(z) 

+ ■■■ + (-T_ lU )z m - 3 G™- 3 (z)G T * (z) + z m - 2 G™- 2 (z)G T »(z) 
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m-2 

= Y,(-T-iur- 2 - k z k G k T (z)G T *(z) 

k=0 

= {-T-lU) G T *(z) . 

( — T-iU) i zG T {z) — 1 

This implies that \Re >n - m \ equals the coefficient of z m ~ 2 in the series expan- 
sion of 

(-T_iu) m - 2 G T * (z) (-r_iu) m - 2 (d + bz) 



(r_iu)- 1 zG r (z) + 1 (r_iu)- 1 z(a + bz) + (1 - zf 

(-l) m (r_in) m - 1 (5 + 6z) 



r_i?i + (a — 2t-\u)z + {b + t^\u)z 2 

Since (a — 2r_iu) 2 > 4t_i«(6 + r_iu), the distinct roots of the quadratic 
equation 

T_iii + (a — 2r_iu)z + [b + r_iu)z 2 = 

are «i and a 2 - Hence, the right-hand side of (9) can be written as 

(-l) m (r^in) m - 1 (q + 62) 
(b + t-iu) {ax - z) (a 2 - z) 

(-l) m (r„in) m - 1 (a + 6ai) (-l) m (r- 1 u) m - 1 {a + ba 2 ) 



{b + T-\u)(ai - z)(a 2 - ai) {b + T-iu){a\ - a 2 )(a 2 — z) 

_ (-^(r-mr-^a + bai) ~ / z \ fc 
(6 + r_iii)ai(a 2 - ai) j^Vai/ 

| (-l) m (r^ 1 u) m - 1 (d + ba 2 ) z \ k 

(6 + r_i«)(ai - a 2 )a 2 ^ \Oi 2 ) 

Thus, we conclude that, for 2 < m < n, 

= (-l) m (T„m) m - 1 (a + fea 1 ) + (-l) m (r-m) m - 1 (a + 6a 2 ) 
(6 + r_iu)a™~ 1 (a 2 — ai) (6 + r_iu)(ai — a 2 )a™ _1 

_ + r_m) m - 2 [a(a^- 1 - a^ 1 ) + 6a 1 a 2 (a^- 2 - <~ 2 ) 

a 2 — ai 

since aia 2 = r_iu/(6 + r_iu). This proves Theorem 1. □ 



3. Inverse. This section evaluates the inverse of Rg tn for some positive 
constant 9. For 1 < i,j < n, if the ith row and jth column of Rg tU are 
deleted, the resulting (n — 1) x (n — 1) matrix is denoted by Rg >n ._ it _j. Then 
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(-lY +j \Re. n;—i,—j\ is the cofactor of the (i,j)th element of the square ma- 
trix Re, n - It is well known (see, e.g., [2], page 582) that 

(10) {Re^= { ~ l)l+ T e T 1 '^ Vl<z,j<n. 

Hence, it follows from Theorem 1 that it suffices to obtain a closed-form ex- 
pression for \Rg n ._i _,-|. We shall, without loss of generality, assume for the 
rest of this section that i < j and i > n—j + 1 since (Re, n )i,j = (Re,n)n-j+i,n-i+i 
and, hence, via symmetry, {Rgl^ij = (i?^) n -j+i,n-i+i- It is convenient to 
define 

2t_iu, if k = -2, 

(1 + w)u - 3(1 + w)u 3 + 2(1 + 2w)u 5 , if k = -1, 
(1 + kw)u k - 3(1 + (k + 2)w)u k+4 

+ 2(l + (k + 3)w)u k+6 , iffc>0, 

r_iu/2, iik = -2, 

((2w- l)u + (l + 2w)u 5 )/2, ifjfe = -l, 
^ (2(1 + kw)u k - 3(1 + (jfe + l)w)u k + 2 

+ {l + (k + 3)w)u k+e )/2, if > 0, 

and for 1 < m < n — 1 , the m x m matrix 

Re,n;-i,-j;m = \\Re,n;-i-j)k,l)\<k,l< m - 

Clearly, Rg jn .-i-j. n -i = jRg jn ._j ) _j. 

Case 1. Suppose i = 1. Using elementary row operations, we observe 
that, for each 1 < m < n— 1, Ro tn --i-j, m can be reduced to the almost upper 
triangular m x m matrix _, ;m , where (-Afl n-— 1,— 7';m)jfe J = whenever 

k>l + 2 and 

(-^0,n;— 1,— i;m)fe,i 

(1 + |Z - fc - l\w)u^ l ~ k -^ , if 1 < k < 2 A m, 1 < / < (j - 1) A m, 
(1 + (/ - k)w)u l ~ k , if 1 < fc < 2 A m,j <l< m, 

r i-k-i, if 3 < k < m, — 1 < I < (j — 1) A m, 
Ti-k, if 3 < k < m, j < I < m. 

Case 2. Suppose i = 2. Using elementary row operations, we observe 
that, for each 1 < m < n — 1, Ro, n --2-j;m can be reduced to the almost upper 
triangular m x m matrix Ag „._2 j _,- ;ni , where (Ag >n .-. 2-j\m)k l = whenever 
k > I + 2 and 

. , /(l + K-lH^- 1 !, ifl<Z<(j-l)Am, 
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i a „ . )ol - f(l + K-3H^- 3 l, if l<K(j-l)Am, 

(a \ _(n-4, if 2<l<(j-l)Am, 

^,«5-2,-iiW3,i-|-_ 3j if j</< m , 

> _f^_jfc_i, if 4<fc<77i,fc-l<Z< (j'-l)Am, 
l^,n;-2 ,-j ; m;*,i - j ^_ fc > ii A < k < m, j < I < m. 

Case 3. Suppose 3 < i < n. Using elementary row operations, we ob- 
serve that, for each 1 < m < n — 1, i?g ,-. m can be reduced to the al- 
most upper triangular mxm matrix Ag jn] _i_j ]m , where (-Afl ) n;-i,-.7^i)fc,j = 
whenever k > I + 2 and 



(^4.^,71;— i,— j ;m) A;, Z 



+ fcHul'" fc l, if l</c<2A?n, 

l<i<(j-l)Am, 
(l + (/-A; + l)u;)u'- fc+1 , if 1 < k < 2 A m, 

j <l <m, 

T;_fc, if 3 < k < (i — 1) A 771, 

(^e,n;-i,-j^i)fc,j = ^ - 1 < I < (j - 1) A m, 

Ti-k+i, if 3 < k < (i — 1) A m, j < I < m, 

n-i-i, if i - 1 < I < (j - 1) Am, 
n-i, if j < I < m, 

fi_i_2, if i < I < (j - 1) A 771, 

fi_i_i, ifj<^<m, 

r;_fc_i, if i + 2 < fc < m, k — 1 < I < (j — 1) A m, 

Ti-h, if z + 2 < A; < m,j < Z < m. 



(^4-0, n;— i,— j';m) i,i 
(^6,n;— i,— j';m)i+l,i 
(^4.0,^;— i,— j;m)fc,Z 



Proposition 2. For 1 < i < j < n,i > n — j + 1, the determinant of 
Re,n;-i,-j;mi 1 < m < n — 1, satisfies the recurrence relation 



\Re,n; 

— ^ ' ( 1) (^4.(9,71;— i,— j;n— l)m— k,rn \Ro,n;— i,— j;m— k~ 1 \ 



-i,—j;m\ 
m—1 



k=0 

fe-1 

x (^.e,n;— i,-j;n— l)m-Z,m-Z-l ; 
Z=0 

w/iere |i?0, n; _; _j ;0 | = 1. 

Proof. Since \Re^ n ;~i-j;k\ = |^.6»,n;-i,-j;fc| for all 1 < k < n — 1, Propo- 
sition 2 follows by using elementary row operations to reduce Ag in ,._j _ J;m 
to a diagonal matrix. □ 
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Theorems 2, 3 and 4 below compute exact expressions for the determi- 
nants \Re^ n --i-j\, where i < j and i>n — j + 1. 



Theorem 2. For 2 < m < n — \, 

\-^9,n\ — 1,— n;m\ 



Proof. We observe from Proposition 2 and Case 1 that, for 2 < m < 
n — 1, 

m— 1 

\Re,n;-l,-n;m\ = ( _ 1)™ ^ 1 (^0,ra;-l,-ri;n-l)fc+l,m|-R0,ri;-l,-n;fc| 
fc=0 

m— fc— 2 

,n\— 1,— n\n— l)m— l,m— I— 1 

— (^4$,n; — 1,— n;n— l)m,m \Rd,n; — 1,— n;m— 1 1 

= |-R0,n;-l,-n;2|T™L 2 

2 2 m-2 
= W U T_\ . 

This proves Theorem 2. □ 

For 1 < m < n — 1, let .Ag n - m denote the m x m matrix such that 



(11) 



( T Tl r 2 

T-l T 



V o oo 



Tm—2 T m — \ \ 
I'm— 3 T m —2 
I'm— 2 T~ m —3 



to n 

T-l T 



J 



Theorem 3. For n — 1 < j < n, 

\R 6 ,n;-2-r,n-l\ = T^* {?- 2 [(l + 2w) 2 U 4 - 1] +f-iT^}\As,mn-j\, 

where t£ is as in (6) and |^4e,n;ol = 1- 

Proof. We observe from Proposition 2 and Case 2 that for n — 1 < j < 



n, 



■2 -Jill - 1> 



|-Rf,n; 

|^,n;-2,-i;2| = (1 + w)u - (1 +ttf)(l + 2«;)u 3 , 
2 

|-Re,n;-2,-j;3| = fc ( J 46»,n;-2,-i;n~l)fc+l,3l^6»,n;-2,-j;fc| 
fc=0 
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,n;— 2,— j;n— 1)3— 1,2— I 

1=0 

= (l + 2w)Vf_ 2 -f_ 2 

+ f_i[(l + u>)u- {l + w){l + 2w)u 3 }, 

\Rd,n;-2,-j;n-l\ = \Rd,n;-2,-j;3\ T -l 4 \Ao,n;n-j\- 

This proves Theorem 3. □ 

Theorem 4. With the notation of Appendix B, for 3 < i < n — 1, 

\Re,n;-i,-i,m\ = \Re,n;m\ VI < 771 < i — 1, 

and 

C 71 (oai + u 2 6) i_ 3 
|i?0 )n; _j_j;j| = — '■ [-{b + T_iu)ai\ l 



(12) ai ~ a2 

-(b + T- 1 u)a 2 ] i ~ 3 . 



CV i2 (aa 2 + u 2 b) 



ct\ — a 2 

For 4<z + l<m<n— 1 we /iaue 

i Ci i ir„iit 4 (aai + u 2 6) 



(13) 



(ai - a 2 ) 2 
x {(ax - n 2 )[-(6 + r_i-u)ai] m - 4 

- (a 2 - u 2 )[-(6 + r_m)ai] 4 - 3 [-(6 + r_i«)a 2 ] m ^ 1 } 
Ci j2 r_iu 4 (aa 2 + u 2 o) 
(ai - a 2 ) 2 

x {(ai -n 2 )[-(6 + r„ 1 n)ai] m ^" 1 [-(6 + r_m)a 2 ] i - 3 
-(a 2 -u 2 )[-(b + T_ lU )a 2 } m - 4 } 

C 2i ir_iu 6 (aai + u 2 6) 



(ai - a 2 ) 2 
x{[-(6 + r_ 1 «)a 1 r- 4 

- [-(6 + r_iu)ai]*- 3 [-(6 + t^u)^]™-*" 1 } 
C 2i2 T_iu 6 (aa 2 + u 2 6) 
(a-i - a 2 ) 2 

x{[-(6 + r_ 1 «)a 1 r- i - 1 [-(6 + r_ 1 «)a 2 r 3 

_[_( fe + T _ lW ) a2 ]m-4 } _ 

vlfeo, /or 3 < i < n — 1, 

|-Re jfl ;-i ,-{i+l);m\ = \Ro,n;m\ V 1 < 771 < i - 1, 

C 3j i(aai + u 2 6) j_ 3 
(14) l-"e,n,-i,-(i+i);il = — _ Q2 [-(6 + r_i«)aiJ 

g3,2(a«2 + M 2 ^) r . v u_ 3 

! [-(6 + r„iu)a 2 J 

ol\ — a 2 
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and for 4<i + l<m<n — 1, we have 

\R$,n;-i,-(i+iy,m\ 

C , 4 i ir_iu 4 (aai + u 2 b) 

(ati - a 2 ) 2 
x {(ai-u 2 )[-(b + T- 1 u)a 1 ] m - 4 

- (a 2 - u 2 )[-(b + T_ lU ) ai Y~ 3 [-(b + r^ lU )a 2 } m - 1 - 1 } 
C^ 2 T-\u A (aa 2 + u 2 b) 
(ai - a 2 ) 2 

x {(a x - u 2 )[-{b + r^u)^]™- 4 - 1 [-(b + r_i«)a 2 ] i_3 
(15) -(a2-u 2 )[-(b + T- lU )a 2 ] m - 4 } 

C^iT-iU (aai +u 2 b) 
(a-i - a 2 ) 2 

x{[-(b + T^u) ai ] m - 4 

- [-{b + T^u) ai } l - 3 [-(b + T- 1 u)a 2 ] m - i - 1 } 
C5 t2 T-\u 6 (aa 2 + u 2 b) 
(ai - a 2 ) 2 

x {[-(& + r„ 1 n)a 1 ] m - l - 1 [-(6 + r^aa]' -3 

-[-(6 + r_m)a 2 r- 4 }. 

Finally, for 3 < i < n — 2, 

C 6 ,irV" 2 (aai + ^ 2 ft) r , w_ 3 

-K0,n;-i,-n;n-l = + r -l u ) Q l 

ai — a 2 

(6 + T_iii)ai2f 



Ce^i* 2 (aa 2 + u 2 6) r ^ i _ A _ i-3 



ai — a 2 
and for 5<i + 2<j<n — 1, 

|^,„;-i,-j ; „-i| - ( ai _ a2 )2 

x {(aai + u 2 6)[-(6 + T. 1 u)a 1 } n - j+i - A 

- (aa 2 + u 2 b)[-(b + r_ 1 u)a 1 ] J - 3 [-(6 + r-jtijaj"" 3 '" 1 } 

C^ 2 t^S\~ 2 (aa 2 + ti 2 6) 
(ai - a 2 ) 2 

x {(aai +u 2 6)[-(6 + r_iu)ai] n ~ i ~ 1 

x [-(b + r^ 1 u)a 2 ] i ~ 3 

- (aa 2 + n 2 6)[-(5 + T- lU )a 2 ] n - j+i ~ 4 }. 

Proof. The first equality is immediate from the definition of R$ in --i —%- m - 
To prove (12), we observe from Proposition 2, Case 3 and (33) that, for 



(16) 
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3 < i < n — 1, 



\Re,n; 



i-1 



,n; — i,— i;n— 1 



-i\k I 



fc=0 



i-fc-2 

X \\ (Ae,n;-i,-i;n-l)i-l,i-l-l 
1=0 

\i— 2 if i— 3- 



(-1J T i-l r -l T -2 
j-2 



fc=2 

\i— 2_* 3- 



(-1)^T*-1^-1 T-2 + r |iV;i-l| 

(-l) i - 1 ril 3 f_ 2 ^- 1 (l - n 2 ) 2 (aai + u 2 6) 



+ 



(ai - a2)[ai(b + r_iii) - v-iii] 
ai(b + r-iu) li 



1 



(-l)^ 1 rr i 3 f„ 2 n'" 1 (l-n 2 ) 2 (aa 2 + ^ 2 6) 
(ai - a 2 )[a 2 (& + r_iii) — r_i«] 

a 2 (6 + r_iu) ri 



1 



T_iU 

r/~ i „,2 



+ 



(-1)«-Vl^f_ 2 [(aai + u 2 ^,! - (aa 2 + n 2 6)5 M , 2 ] 



ai — a 2 



where Si m j is as in (31) and (12) follows after some algebra. For 4 < i + 1 < 
m < n — 1, 



m— 1 

— i,— i;ml ^ ' ( 1) ~ (^0,n; — i,— i;n— l)fci+l,m \R-9,n; — i,— i;k\ I 
fci=0 

m— fci —2 

x (^4.0,n;— i,-i;n-l)m-J,m— /— 1 



=0 



m—1 



+ (-T-l)" 1 1 X! r m-fci-l(- r -l) kl \R6,n;-i-i;kl\i 
k\=i+\ 
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where 



i-2 



(17) 



im = T4_iT_i'f_ 2 f_i + ^(-1) ^ 1 T fn „ fc T_ 1 fc 3 f_ 2 r_l|-R6»,n;A 

+ (-l)~ 4 f m _jT~{'~ f_l|i?g )n;i _i| 

+ ( — 1) * Tm-i-iT_l \Re )n --i-i\i\ 

= T* l _ 1 rrif_ 2 f_i + (— \)~ % T m -iTZ\~ f_l|-Rfl )n; i_l| 

+ (-1) * T m -i-\T_{ \Rg^ n --i-i-i\ 
.ni-l/l „,2\2/'~ i „,2 



+ 



r£ 1 (ai - a 2 )[ai(6 + r_iu) — T_ru] 
f-af-m^H 1 ~ n 2 ) 2 (Q«2 +u 2 b) 
Ti 1 (ai - a2)[a 2 (& + T_iit) - T-iu] 



a±(b + r-i-u) 

r_iu 
a 2 {b + r-i-u) 
r_iu 



i-3 



i-3 



_ f_ 2 f_i[(aai + u 2 b)S i<mt \ - (aa 2 + « fe^.m^] 
r£ 1 (o:i - a 2 ) 

The last equality uses (33). An immediate consequence is |i2g n; _i ) _i.j + i| = 
(— r_i) l+1 £j_|_i. Then repeating the above argument, we obtain, for 4 < i + 
1 < m < n — 1, 



\Re,n; 



(-r-i) m e, 



m— 1 



+ (-r_ 1 r- 1 e ^^(-t^)-* 1 

fcl=i+l 



(-^l) fcl £ fcl 



fel-1 



+ (-r_i) fcl 1 ^ r fcl _ fe2 _i(-r_i) fca |Jfy jn ._ <> _ i . fca | 

fc 2 =i+l 

m— 1 

(-r_i) m ^ m + (-r_i) m_1 X] 

fci=i+l 

m— 1 fci— 1 
+ (-r_ 1 ) m " 2 E E 

I'm— k\ — \T~k\— ki — 1 

fcl=j+2fc 2 =i+l 

X (-T-iJ-^li^^-i-ijfcal 
m— 1 

(-r_i) m £ m + (-r_i) m_1 E T m-kx~\t,k\ 

k%=i+l 

m—l k\—l 

+ (-T-l) m ~ 2 E E T m-fci-l r fci-fc2-l£fc 2 
fci=i+2fc 2 =i+l 
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m— 1 fci — 1 k2 — 1 

+ (-T_i) m " 3 ^ ^ ^ 7-m-fei-lTfei-fe 2 -l 
fcl=i+3 fc2=i+2 fc3=i+l 

x 'Tfe — fc3-lCfc3 + ' " ' 

+ (-r_ 1 r+ 2 

m— 1 fci— 1 fem-i-3— 1 

X "' Tm-fci-lTfci-fca-l • • ■ 

k i =7tl— 2 k2=m—3 fc m _i_2=£+l 



'0 

Hence, it follows from Lemma A. 4 (see Appendix A) that |-Rfl n ;-i,-i;m|)4 < 
i + l<m<n — 1, equals the coefficient of z m in the series expansion of 

(-r_iu) m G{(z) + (-r„ 1 n) m - 1 zG r (z)G 5 (z) 

+ (-r_iu) m - 2 z 2 G2(z)G € (z) + • • • + (-r-iu)^ 1 ^-^^- 1 ^)^^^ 

m— t— 1 

= X (-r-iu)™- ***G*(*)G C (*) 
fc=0 

m— i— 1 

= (-T-iu) m G t (*) E (-r-xuJ-^GjCz) 

fc=0 

= (-r_iu) &£(z) — . 

(— r_iuj i zG T [z) — 1 

It follows from Lemmas A.l and A. 4 that \Rq re ._j _j. m |,4 <i + l<m<ra — 1, 
equals the coefficient of z m in the series expansion of 

(-r_m) m G ? (z) _ ^(-r-iurCoj + b^) 



r_iu)- 1 zG r (z) + 1 (r_iu)- 1 z(a + fcs) + (1 - z) 2 

(-l) m ^ +1 (r_m) m+1 (ag + bgg) 



r_in + (a — 2t_iu)z + (6 + t_\u)z 2 



Letting a\ and a2 be as in Theorem 1, the right-hand side of (18) can be 
written as 

(_l)m 2 i+l( r _ lU )m+l( a? + 
(b + T-iu)(ai - z)(a 2 - z) 

_ (-l) m (r_itt) m+1 (a g + jgai) ^ z fc+ * +1 
(6 + r_i'u)ai(a2 — ai) ^ 
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(-l) m (r-m) m+1 (ag + ka 2 ) ^ 
(6 + T_iu)(ai - a 2 )a 2 ^ a\ 



Thus, we conclude that, for 4<i + l<m<n — 1, 
_ (-l) m (r_m) m+1 (a c + 6 5 ai) 



(6 + r_iu)a™ J (a 2 - «i) 
(-ir(r-mr +1 (ag + ^a 2 ) 

(6 + r_in) (ai — c^a™ -4 
(-l) m (r„iu) i+1 o c 



+ 



(19) = ^ ^'~ ia > - n 2 )[(6 + r^u) ai } 



m—i—l 



&l — OL 2 

-(a 2 -u 2 )[(b + T- 1 u)a 2 } m ~ i - 1 } 
| (-l)™(T-m) i+ V(ag + &g) 
cti — Oi 2 

x {[(ft + r.x^ai]" 1 -^ 1 - [(ft + r.m)^]" 1 -^ 1 }- 



Equation (13) follows from (19) and Lemma A. 4 after some algebra. Next, 
from the definition of ify / i+ ]V m we have 

\Re,n;-i,-(i+l);m\ = \Re,n;m\ V 1 < m < Z - 1, 



and 

\Re,n;-i,-(i+l);i 



i-l 

X] 1 (Ae,n;-i,-(i+l);n-l)k+l,i 

k=Q 

i-k-2 

x \Re,n;-i-(i+l);k\ Yi 
(=0 

(_ + (j _ i^uVfr^ 

+ (-ir 2 [l + (i - 2)zt;]^- 2 rr i 3 r_ 2 
i-2 

,n;— i,— (i+l);fc| 

fc=2 

,n;— t, — — 1 1 

(-irV-^iff-a + f_i|i2fl, Bii _i| 

+ l ^(-iy- k -\_ k _ 1 r!_- 1 k - 2 f. 2 \R e , n . k \ 

k=2 



GAUSSIAN RANDOM FIELDS 



17 



+ 



+ 



(-l)*-ViT- 2 ?/- 2 (l - u 2 ) 2 {a ai + u 2 b) I r (b + r_m) ai 
r_i(ai — 02) [(6 + r_i«)o!i — r_iit] \[ t_\u 

(-l) i - 1 Til 2 f_2« < - 2 (l - u 2 ) 2 (5a 2 + n 2 6) 
(-l)*-i T i"i 2 f_ 2 [(aai + « 2 6)5 M _i,i - (aa 2 + « 2 6)5 M _i ;2 ] 



1 i-3 



(6 + r_iit)o:2 



i-3 



ai — «2 

The last equality uses (33) and (14) results after some algebra. For 4 < 
i + l<m<n — 1 , we have 



m—l 



\ R e,n;-i-(i+l);m\ — ^2(-^) m k 1 (^0,n; 



'i,-(i+l);n-l)k+l, 



k=0 



m— fc— 2 

x \Re,n;-i,-(i+l);k\ Yl (^0,n;-i,-(i+l);n-l) m -l, m -l-l 
1=0 

(-l)™-\l + W )(l +mW )u m+1 T™i 4 T-2f-2 

+ (-l) m " 2 [l + (m - lHu m -VY 4 r- 2 f- 



fc=2 



2 

m— fc— 3 - 



T m -k\ft6,n;-i ,-(i+l);fcl r -l r -2^ 



1 / i\m-i- 175 |„_m— i— 1~ 

+ 1-1) ^m-i|-K0,n;-i,-(i+l);i-l| T -l r - 



+ (-!)'" J f r7 i_j_i|i?0 jn ._ i _(j +1 ). 



„m— i— 1 



il'-l 



m—l 



+ E (-1) 
fc=t+i 

(-r-i) m |„ 



m— A;— 1 



7m-fc-l|-R0,n;-i,-(i+l);fcl r -l 



m— A;— 1 



m—l 



+ (~ T -l) m 1 E T m-fc-l(-T-l) k \Re,n;-i,-(i+l);k\, 
k=i+l 



where 



Cm. = 7",n-l r -l^-27 : -2 + (-1) ^m-il-R^nji-l k-1 ^-2 



-i-1 



il'-l 



i-2 



+ E( _1 ) fc lr m-fc|-Re,n;fck_l fc 3 T- 2 f_ 2 
fc=2 



T_ 2 T-2 + (-l) l fm-i\Re ,n;i— lr— 1 T - 
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+ (— 1) 1 1 f m _i_i|-Re )n; _j i _(j +1 ) ; 



-i-1 



il'-l 



+ 




i-3 



i-3 



— a2)[ai(6 + r_iu) - T-iu] 

f^ 2 f ^U^ 1 {I -U 2 ) 2 (ha 2 +U 2 b) 
_ Q!2)[Q!2(& + T_l«) -r_i«] 

_ T„ 2 T-2[(aQi + u 2 b)Sj, m ,i - (aa2 + tt 2 b)5j, m ,2] 
r^^ai - a 2 ) 

The last equality follows from (33). An immediate consequence is that 

\ R e,n--i-(i+l);i+l\ = (— 

Now repeating the above argument, we obtain, for 4<i-\-l<m<n — 1, 

m— 1 

l^6»,n;-i ,-(i+l);ml = (- T -l) m ^m + (-T-l)"'" X! T m-fci-llfci 

m— 1 fci— 1 

+ (-T_i) m ~ 2 ^ ^ T m-ki-lTki-k 2 -\t,k 2 
fci=i+2 fc 2 =i+l 



+ ... + (_ r _ 1 )*+l T m-i-l e 



Hence, it follows from Lemma A. 4 that l-R^n;— i,-(i+i),ml?4 < i + 1 < m < 
n — 1, equals the coefficient of z m in the series expansion of 

(-r_m) m G f - 



(r_ 1 u)- 1 zG r (z) + l' 
Following the proof of (19), we conclude that, for 4<i + l<m<n — 1, 

|-^0,n;-i,-(i+l),ml 

(-l) m (6 + r_iu) m u 2 ( <+1 > Maf"* - aS 1 -*) + u^a™^ 1 - a^"* -1 )] 



ai — a 2 



«1 — «2 
2> 



x {(ai — u )[(6 + r_iit)ai] 



m— i— 1 



(a2-u z )[(b + T-iu)a 2 } m - i - 1 } 



+ 



(-l) m (r_m) l+1 u 2 (a f - + 6 f -) 



ai — a 2 



x {[(6 + r_i«)ai]'"- < - 1 - l(b + r_ 1 u)a 2 ] ro ^- 1 } 



m— i— 1 ■ 
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and (15) follows from Lemma A. 4 after some algebra. Finally, we observe 
from Case 3 that, for 5<i + 2<j<n, 

|p I _ -J-i-2| R 1 1 A I 

\ r, '9,n;—i,—j;n—l\ '—1 \ rt '9,n;— i,~ j;i+l\\ y ~ L 9,n;n— j\ ■> 

and (16) follows from Lemmas A. 2 and A. 3 (see Appendix A). This proves 
Theorem 4. □ 



Remark. We wish to add that the exact results of Sections 2 and 3 
have been checked by Mathematica [15]. 

4. Asymptotic approximations. This section establishes asymptotic ap- 
proximations for \Re >n \ and Rg^- These approximations are very sharp in 

that the error is of the order 0((2 + v / 3)~ n ) as n — > oo. The following the- 
orem gives the approximation for the determinant. 



Theorem 5. Let < 0o < Pi < oo. Then 



l-Re.nl = [-(6 + r_m)ai] n 2 



a\ — a<i 



w 



3n 



-4 e -2(n-2)»^3 



1 + 
2(2 + v / 3) 1 ™" 1 



1 + 



2- V3 \ n 
2 + v^ 



(12 + 7V3) 
60 + 30^ 



+ 0{w 2 ) 



as n — > oo uniformly over 8 € [/3q, (3\ 



Proof. From the definitions of ot\ and a 2 in Theorem 1 and using the 
Mathematica command 



we observe that 
(21) 



Simplif y[Series[a!/Q; 2 , {w, 0, 0}]], 
a 2 2- v / 3 r 



ai 2 + y/3 



[1 + 0(0], 



as n — > oo uniformly over G [/3o , /?i ] • Using Mathematica in a similar way, 
we also have 



aa<i 



+ u 2 b 2-V3 



l + 0(w)], 



act\ + u 2 b 2 + s/3 
as n — > oo uniformly over £ [/3q, Thus, it follows from Theorem 1 that 



aai + u 2 b 2 
-Rfl.n = [-(6 + r_iu)ai] n 

«1 — «2 



l + O 



2--y/3 

2 + v^ 



20 



W.-L. LOH 



as n — > oo uniformly over 6 £ [/?o,/3i]. Next, using the notation of Theorem 
1, we observe that 



e 4w (a-2r_m) = 1 + 4we 



2w _ g4ui 



8w 3 ~ 3!2 fc (2 fe+2 -A;-3) 



E 



w 



k=0 



and 



e 8w [(a - 2t_ 1 u) 1 - 4r_iu(6 + r_i«)] 



(A + 3)! 



(1-C« 



where 



c« = 96 ™ E ti^t - 18432w E 77^7 + 209952w E 



» (4w) fc ~ (Q w )k 

1536 ™Et7^ + 5832w E 



fc=0 



= -4tu(l + 0(w)), 
as n — ► 00 uniformly over 9 £ [/?o,/3i]. Hence, 
(6 + r_i«)ai 



(22) 



2 

-4ui 



[e to (a - 2r_i«) + e 4 ™ J (a - 2r„ lU ) 2 - 4r_ lU (6 + r_i«) ] 



4w 3 ^ 3!2 fc (2 fc + 2 - k-3)w k 2w 3 1/2 
~3~ ^ (* + 3)! + 73 ( " Cu ° 



k=0 



-Aw 



4w 3 3!2 fc (2 fc+2 -fc-3) 
2io 3 



10" 



+ 



V3 



(fc + 3)! 



' j ^ 4 *+ifc!(fc + 2)!_ 



2(2 + a/3)io 3 



(1 + Ai,, 



where 

Ai (l0 



2^-^ ~ 3!2 fe+1 (2 fe+3 -£;-4W fc A 4? „ ^ (4w) fc 



2 + V3,e 



(fc + 4)! 



-2w + 
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60 + 30^ 



+ 0(w*), 



as n — > oo uniformly over 9 G [/?o,/3i]. We conclude from (22) and Theorem 
1 that 



(2 + V3) 



w 



3n-4 r 



2(2 + V / 3)(l + Ai,, 



2^/3 

2(3 + V3)w 
3(2 + ^3) 

w 3n-4 e -2(n-2)«,^3 



n-2 



+ 0(w' 
2(2 + ^) 



n-1 



1 + 



(12 + 7^3) nw 



+ 0(V 



60 + 30\/3 

as n — > oo uniformly over £ [/3o,/3i]. This proves Theorem 5. □ 

The next two theorems give the approximation for the inverse Rq\- 

Theorem 6. Let < (5q < (3± < oo. Then with the notation of Appendix B, 
sn-3 to{t- 2 [(1 + 2w) 2 u 4 - 1] + f-irfKai - a 2 ) 



(Rgl)2,n-l - (-1)' 



u 
a i 



n-2 



l + O 



T^ 1 (aai + u 2 6) 

( 2-Vs y 

V2 + v/3/ 



(^n) 



n—l,n—l 



C 7 ^(aa 2 + u 2 b) 



a 2 



(6 + r_iu)ai] 2 + u 2 6)[-(6 + r_iu)ai] 2 \ai 

l + O 



n— 4 



(^n)En = (-l) 



( 2- V3 \" 
V2 + ^3/ 



w+1 -w 2 ?x 2 (ai - a 2 ) 



n-2 r 



l + O 



(02,n = (-l) 



n+2 



r_i(aai + u 2 b) V «i 
{f_ 2 [(l + 2u;) 2 n 4 - 1] + f-i^Xai - a 2 ) 



/ 2-yg y 
V2 + v/3/ 



r 2 1 (aai + u 2 6) 



n-2 



l + O 



(( 



2- 



,2 + ^3 

(-lf+ic 6il r^r- 2 [-(& + r_iu)ai]- 3 



[-(6 + r_iu)ai] 



n-2 



l + O 



/2-V3V 
V2 + x/37 
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jn-i— 2(%„,_ i „,2l 



( R 9,n 



(-lf^Ce^ (aa 2 + u 2 b)[-(b + r-m)a 2 ] 
(aai +u 2 6)[-(6 + T_iu)ai] n - 2 

1+0 "»-vsv* 

C3,l 



i-3 



-(b + T-iu)ai] 



2 + ^/3 
l + O 



1 



i,i 



(b + T-iu)a\ 
1 

(b + T_iu)ai 
C31 



l + O 



( 2-VS 
2-V3\ n 



2,1 



-(6 + T-iu)a\\' 

i-3 



l + O 



2 + V3 



(- C 6>1 r^ [- (6 + r_i«)ai]"- 



■i-2 



l + O 



[-(6 + r_iu)ai] n - 2 
(-l) l+1 06,2^(002 + u 2 6)[-(6 + r-iujaa]"-*- 2 



/2-V3V 
V2 + V3/ 



l + O 



(aai + n 2 6)[— {b + r_iu)o!i] n " 
V2 + y/3j 



and 

( R 9, n )n-l,l = (~l) n+2 



{f_ 2 [(l + 2-u;) 2 n 4 - 1] + f_i7f}(ai - a 2 ) 



t 2 1 (aai + u 2 b) 



u 



n-2 



1 + 



( 2-V3 \ n 
K2 + V3J 



as n — > 00 uniformly over £ [/?q, /3i] and 3 < z < n — 2. 

Proof. Theorem 6 is a consequence of (10), (21) and Theorems 1-4. 



□ 



Theorem 7. Let < (3q < (3\ < 00. Then with the notation of Appendix B, 



a\ — u 



a 2 — u 



a 2 



Cl } lT-lU 4 

ot,\ — a,2 I [— (b + r_iii)ai] 3 [— (6 + r_iu)ai] 3 V«i 

Oi ) 2T_iu 4 (aa2 + u 2 b) 
{a\ — a 2 )(aai + u 2 b) 



n—i—2 
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a± — u 



0: 9 



Cl2 — U 



[-(6 + r_iu)a 2 ] 3 [-(6 + r_iu)a 2 ] 3 \ai 



0'2 



n-2 



+ 



1 



1 



0: 2 



C 2 ,ir_m 6 

— a 2 I [— (ft + r_in)o;i] 3 [— (6 + T_iu)ai] 3 \ai 

C 2t 2T-iu 6 (aa2 + n 2 6) 
(ai — a 2 )(aai + u 2 6) 

1 



a 2 



1 



a 2 



1 + 



[-(6 + r_iw)a 2 ] 3 \ai / [—(6 + T_iu)a 2 ] 3 \ai 

p-\Z3 V 
V2 + v^/ 



as n — > oo uniformly over 9 £ [/3q, /3i] and 3 < i < n — 2, 



" ri^ax - a 2 ) V«i 
x (c 6 ,i 



i-i+2 



(aai + u 2 o) — (aa 2 + n 2 o) ( — 



n-j-1 



C 6j2 (da 2 + u 2 b) 
aai + k 2 o 



2 \ ( a 2 

(aai + n o) — 



i-3 



(aa 2 In 



a 2 



n—j+i—i- 



l + O 



2-\/3\ n 



VV2 + ^/3. 

as n — > oo uniformly over G [/?o, and 5<z + 2< j<n — 1, and 
(-^nki+i 

C4 lT_i1i 4 



ai — a 2 



ai — u 



a 2 — u 



a 2 



[— (6 + r_iu)ai] 3 [— (6 + r__iu)ai] 3 \ai 
C4 )2 r_iu 4 (do; 2 + ii 2 6) 



n— t— 2 



(ai — a 2 )(aai + n 2 o) 
ai — u 2 



a 2 



a 2 — u 



[— (6 + r_iu)a 2 ] 3 [— (6 + r_iu)a 2 ] 3 \ai 



0'2 



n-2 
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+ 



ol\ — a 2 



1 



1 



«2 



-(b + r_iu)ai] 3 [— (6 + r_iii)ai] 3 

C 5i 2T-iu 6 (aa 2 + u 2 b) 
{a\ — a2){aa\ + u 2 b) 



n—i—2 



A 9 



1 



a o 



[-(b + T_ 1 u)a 2 \ 3 \ai/ [-(6 + T_iu)a 2 ] 3 \«i 
2->/3\ n ' 



n-2 



l + O, 

,2 + ^3. 

as n — > oo uniformly over 6 [/3cb /Si] and 3 < i < n — 2. Finally, 



, n—i—2 

-(6 + T_in)ai] 3 \ai) 

C 6i2 (aa 2 + n 2 6) 

(aai + u 2 6)[-(6 + T„iii)Q 2 ] 3 V«l 



u \ n ~ i - 2 (a 2 



a i 



,2 + V3, 

as n — > oo uniformly over £ [/?o, /3i] and 3 < -i < n — 2. 

Proof. As in the proof of Theorem 6, Theorem 7 is a consequence 
of (10), (21) and Theorems 1-4. □ 

5. Fisher information. Stein [12], page 178, noted that the asymptotic 
theory of maximum likelihood estimation based on independently and iden- 
tically distributed observations suggests that calculating the Fisher infor- 
mation matrix is a fruitful way of learning about the behavior of the max- 
imum likelihood estimators. Abt and Welch [1] showed via three examples 
(either proven analytically or justified by simulation) that, for the covari- 
ance parameters of Gaussian processes under fixed-domain asymptotics, the 
covariance matrix of the limiting distribution of their maximum likelihood 
estimators equals the limit of the inverse Fisher information matrix. Also, 
Section 6.6 of [12] reported a numerical study of the Fisher information ma- 
trix for observations from a mean Gaussian process on the real line with 
a Mater n- type covariance function. 

Let X n be the random vector as in Section 1, having a n^-variate normal 



distribution with mean and covariance matrix £„ 



■ n . Motivated by 
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the work described in the previous paragraph, the following theorem com- 
putes the Fisher information matrix for the parameters <f),0i, ■ ■ ■ ,9 d based 
on X n . 



Theorem 8. LetX n ~ JV^fOjE^j] , (J. Then the elements of the Fisher 
information matrix for the parameters (f), #i, . . . , 9 d are given by 



-E 



2 



logL n ((/),9 1 ,...,9 d ) 



-E 



del 



logL n (ct>,9i,...,9 d ) 



-E 



-E 



2 



log L n (<j),6i,...,9 d ) 



d 2 n d 
2^ ' 

n d - 1 {29 t + 5) 
dn d ~\9 t + 2) 



+ 0(n d ~ 2 ), 
+ 0(n d - 2 ) 



2 



■log L n ((j),9 1 ,...,9 a 



de s de t 

as n — > oo whenever 1 < s ^ t < d. 



0(ri 



d-2\ 



PROOF. We observe from [2], page 600, that 



(g)%,n 



t=l 



Y[[\ R e t ,n\ 



t=l 



and, hence, it follows from (5) that 
2logL n (^9 u ...,9 d ) 

^ log(2vr) - dn d log - dn d log(<£) + 3n d ]T ^g(6 t ) 



(23) 



-n 



t=i 
d \ -1 



n d - l Y.\°g{\Re t ,n\) ~ 2de } d ',f d K[ ® ) Xn. 



t=l 



7T" 



t=l 



Consequently, 



-E 



2 



logL n ((f),9i,...,9 d ) 



dn d d{d+l)n d _d 2 n d 
2^2 + 2§ 2 2(j) 2 ' 



Writing wt = 9t/n, we observe that 







3n a 



aSj[3n <, M0t)] = — ar, 



t) 2 



(90 



t 

d-l 







't 



(24) ^(n-logl^l) m 

n d ~ 3 d 2 \R 0un \ 



d\Re un \ 

„d-3 



f d\Re un \ 
\Ro un \ dw 2 \Re un \ 2 \ dw t 
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Next, we observe from Lemma A. 6 that, as n — > oo, 

a_ 



(25) 



1 

-tr 

n 

1 

-tr 

n 



(w d R e d ] n ) R ^ 



12 2{A9 d + 9) 



d 2 



n 



+ 0(n 



We observe that 
2 



01 " • ' #d^n ^ (§) -Rft ,nj A n 



d 



■d-l 



t=l 



ml--e d x' n i^R^ n )x n 



. t=i 



+ 66l-.6jX' n 



+ 0-i ••• o d x' 



/d-l 

® 

\ t=l 

d-l 



. t=i 



X r , 



X r , 



and, hence, using (25) 
d 2 



E 



de 2 



IT 



dAd 



9t ,n 



Xr, 



. t=l 



6(2 



d\Q'i Q 
1 - ' • "d 



7T 



d,Ad 



tr 



(26) 



, 6(2")fl?---flj tr 



0^1i ) S *.9i,-.»<i;n 



(=1 



r dAd 



+ 



7T' 



'd-l 



7T 



d^d 



tr< 



g n d g^d 1 

2n d - 1 (2fl d + 3) 



'd-l 



t=l 



( R eln) 



t=l 



dd 2 



( R eln> 



S 



+ 0(n d ~ 2 ). 



+ n^ 1 tr 



<t>fii,...,O d 
{d0j Ke < 
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We conclude from (23), (24), (26) and Theorem 5 that, as n — > oo, 



d 2 



3n d (3n d 2n a 



d-l\ „d-l 



29 2 \29 2 e 2 



+ 



n d ~ 1 (29 d + 3) 



+ 0(n 



d-2\ 



n d - 1 {26 d + b) 



+ 0{n d - 2 ). 



Finally, using (25) again, we have, as n — > oo, 

d 2 



-E- 



89, 



log L n ((j),9 1 ,...,9 d ) 



-E 







2 d - 1 d9 z ■ ■ -9 3 ~ , ( d \ 1 ~ 

n d0d+l dX 'n ( Rd t,n ] X n 



d9 d 

3(2 d - 1 )d9f---9j 



t=\ 



7T 



dAd+1 



tr 



<g)Ro t 



n z->(t>fii,...,6 d ;n 



2 d ~ 1 d9f---9 d 

it-d^d+i 

3dn d dn d ~ x 



tr< 



2cj>9 d 20 
dn d ~ l {9 d + 2) 



tr 



L \ t=l 
- /d-1 

Ai=l 




Rq „ \R 



+ 0(rt 



d9 d "^ n > 

d-2\ 



and 



89 d _ 1 d9 ( 



log L n ((f),9 l ,. .. ,9 d ) 



E- 







d9 d .! 



3(2 d -i)9f---9 2 ~, 



7T 



lR e t ,n ) X n 



+ E 







d9 d ^ 



t=i 
d-i 



7T 



dMd 



t=l 



9n d ?>n d - 1 

+ -r^— tr 



29 d -\9 d 29 d 



d9 d -! 



+ 



3n 



d-l 



29 



■tr 



d-i 



gg R e d )n) R O d ,n 
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n d ~ 2 \( d 



R e d \ 1 ,n) R dd-i,n 



tr 



0(n 



d-2\ 



The proof of Theorem 8 is completed by invoking the symmetry of 9\, . . . , 9 d . 

□ 

6. Estimating the scale parameter. Let X n ~ N n d(0, q 1 ...,9 d -n) be as 
in Section 1. Also let < /?o,t < fii,t < oo,i = 0, . . . , d, be known constants 
such that /3o,o < <fr < /?i,o and (3o,t < 9t < 0i,t,t = l,...,d. This section will 
be concerned with the construction of a consistent maximum likelihood-type 
estimator for the scale parameter (j). First set the estimator for (6\, . . . ,64) 
to be a known arbitrary but fixed vector, say (9i,...,9~d) £ Ilt=i(A),i)/0i,t), 



and define the estimator 



h,...,6 d 



for d> to be that value of 6 that maximizes 



the log-likelihood function log L n (<p,9i, ... ,9d) [see (23)]. On differentiating 
the log-likelihood with respect to <j>, we obtain 



d 



dn a 



— [2 logL n (<f>,0 1 ,...,e d )] = - + 



d2*9l--9 d ~, 



7T 



dAd+l 



x„ 



9t,n 



x n . 



. t=l 



Equating the right-hand side of the above equation to zero and then solving 
for (j) gives us a maximum likelihood-type estimator for <f>, namely 



2 d el--9 d ~, 



TT d n d 



-1 i/d 



t=l 



Theorem 9. With the above notation and conditions, we have 



2 d 9f' 



■01 



TT d n d 



-E 



\ t=i / 

d \ 



Xn 



+ 0(n- 1 ), 



as n — > 00. Consequently, 0^ ^ — > </> in probability asn^oo. 

Proof. Since (#1, . . . ,9^) is a constant vector, we observe from [2], page 
600, and Lemma A. 5 (see Appendix A) that 



2 d 9j • • • 9 6 d 
n d n d 



E 



x' 



Xn 



t=l 
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<P d 6f ■ ■■e d 
ef--- e 3 d n d 



tr 



■ d \ -1 / d > 



L \ t=l 
1 



. t=l 



tr^r 1 R eun ) 



= ^ + 0(n- 1 ), 
as n -> oo. Next, we observe from [8], page 53, that 



2 2d e\ ■ ■ ■ e d 



£7 



i? 



6»t,n 



2^> M #! j ) 



ef---e d n 2d 



tr 



2n 



e t ,n 



k ,n 



+ 



M 0f ■ • • flg 



0? • • • e d n 2d 

Hence, it follows from Lemma A. 7 that 



. t=l 



Var< 



t=i 



(27) 



tr 



TT d n d n I 

_ 2cj) 2d ef ■ ■ ■ e d 
~ • • ■ eW d 



0(n- d ), 



R e!,n Re ^ 



as ?i — > oo uniformly over /3o,o < 4>< /3i,o and /3o,t <0t,0t < Pi,t,t = l,...,d. 
The final statement of Theorem 9 follows from the definition of 0^ g , 
Chebyshev's inequality and the fact that the dth root function is continuous. 

□ 



Remark. Theorem 8 shows that the amount of Fisher information con- 
tained in the sample on the scale parameter <f> is an order of magnitude 
greater than that on the correlation parameters 0\,...,0d- Thus, it should 
not be really surprising that incorrect specification of the values of 0±, . . . , Od 
can still lead to consistent esimation of cj>. Crowder [4] (page 49) has a dis- 
cussion of such a phenomenon in a more general setting (see also [12], page 
175). 
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7. Sieve maximum likelihood estimation. Let X n ~ N n d(0, S ( ^ i g 1: ... ) g d;n ) 

be as in Section 1. The following definition is taken from [12], page 163. 

Definition. For a class of probability models {Pg for a random 

field on a given bounded domain, a function h on S is said to be microergodic 
if and only if, for all £, £' 6 S, 7^ /i(£') implies Pg is orthogonal to . 

We observe from [12], page 164, that, for d=l, <f> is microergodic, while 
9\ is not. This is consistent with Theorem 8 which states that the Fisher 
information for Q\ is bounded as n — > 00 if d = 1. We shall assume in this 
section that d > 3. 

Let < v < (d — 2)/(d + l) and < /? ,t < /?i,t < 00, £ = 0, . . . , d, be known 
constants such that /3q,o < 4>< Po,i and /?o,t <9t < 0i,t, t = 1, . . . ,d. We define 
a sieve f2 n on the parameter space of (0, 9±, . . . , 9 d ), namely, 

n n = \ ( ~) • ■ • ; —7. ) -l3o,t<^: < Pi,t,k integer, < t < d\. 

The sieve maximum likelihood estimator for ((f), 9\, . . . , 9 d ) is that element 
(0,01, ... , 9d) G such that 

L n (4>, 0i, . . . , d ) = sup{L n (</>, 9x,...,9 d ):(4>,9x,...,9 d )e O n }, 

where the likelihood function L n is as in (5). For each e > sufficiently 
small, define 

fi nj£ = O n \0 - e, + e) x (0i - e, 9i + e) x ■ • • x (9 d - e, 9 d + e). 

Theorem 10. Lei d> 3, X n &e as in Section 1 wit/i covariance matrix 
^(j>,e 1 ,...,e d ;n- Let (<p,9i, . . . ,9 d ) be a constant vector (depending only on n) 
such that (cp, 0i, ... , 9 d ) £ Q n , \<p — <f>\< n~ u and \9t — 9t\ < n~ v , t = 1, . . . , d. 
Then for any e > 0, 

lim p r su P {L n (^, § u ...,e d ): {4>, 0i_, ...,e d )£ n w , e } <e \ = 1 

where L n is the likelihood function given as in (5) . Consequently, the sieve 
maximum likelihood estimate ((f), 0i, . . . ,9 d ) — > (0, 0i, . . . , 0^) in probability as 
n — ► 00. 

Proof. We observe from (23) that 

2 - - - 

—j [log (0, 0i , • • • , 9 d ) - log L n (</>, 0i , . . . , 9 d )] 

df---o 3 d 
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0t,n I ^™ 



ir a (b a n a 



6t,n 



X n . 



t=i 



Now, 



P\ sup {log[A»(^ei,...,^)]-log[L n (^Oi,...,e d )]}>log(e; 

l,e 1 ,...,e d )en n , e 



pj |J {log[L n (<Mi,---,#d 



]og[L n $,5 1 ,...,e d )]}>]og(e) 



< p{iog[L n {4>,h,---M] 



(4>,e 1 ,...,e d )en nt£ 



-io g [L n (4>,e 1 ,...,e d )]<-iog(e)} 

\ n d 4> d n d n \^[ 



R 8 t ,n X 'n 



t=l 



E 



od/T3 23 / d 



n dAd n o 



x n 



t=l 



— d -r d - d -x n Q9 u eun x n 



\t=i 

<-log( 4- I - log I J" J 



--.L lo r 



\Re 



t=l M"0 t ,n 

2 d e\ ■•■el 



TT d (b d n d 



EX' n \®R~ et \ X n 



•32 



W.-L. LOH 



+ 2de !::- e / Ex'J <^R f 



r K d (b d n d 



21og(e) 



71 



Hence, using Chebyshev's inequality, we obtain 



P{ sup {log[L ft (^e 1 ,...,e d )]-log[L n (^ > ei,...,5 d )]}>log(e; 

'Opi,...,0 d )en„, e 



< 51 Var< 

(^,0i,...,e d )eOn,s 



TT d (h d n d 



t=i 



X r , 



2d9f "' e3d -x'J 6dm 



TT d d> d n d n 



,n ^« 



(28) 



x | - loi 
1 



3 log 

'\ R eJ 



t=i / 



n t=l \\ R 6t,n\ 



odfl3 a3 / d 

+ 2& K:' e < EX'J 6dR 



TT d (b d n d 



e t ,n X n 



21og(e) 



t=i 



Next, we observe from Theorem 5 and Lemma A. 5 (see Appendix A) that 



log ^ +31o 



n ■ ■ ■ vd 



71 ^— ' 



.( \ R ~eJ 



+ 



2 d ef ■■■el 



EXL 6dRs J X n -1 



n d 4> d n d n VSi 



3n - 4 . J 



'i * * * ' 



'i ■ • ■ o d 



t=i 



0A 2(n-2) 



+ ;rMr 



<g)(R~ t ] n R 9t ,n) 



t=l 



l + 0(n~ 2 ) 
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t=l 



t=l 



+ 0(n~ 2 ) 



loe 



i=l 



n ~ \ Ot I n f— f 



n - vt 



t=i 



t=i 



^ I +- 

n 



3&/0A 4 



4 U 



+ 



t / p t 



+ 



+ 0(rT 2 ) 

Ad \ Ad d 



as n — ► oo uniformly over /?o,o < < A),i and /?o,i < t < j3\,u t = 1, . . . ,d, 
where 

/(i)=t-log(t)-l Vf>0, 



<7(0t,0 t ) = -41og 



2(0* - t ) + 



+ 



3^ 6t{0 t \ dtfO t 

We further observe that f(t) > for all t > 0, /(i) = if and only if i = 1, 
4r9WtA)>0 V6 t >e t , 

4r9(PtA)<0 ve t <e u 
o6 t 

g(9 t A) = if and only if t = 0+ 
Hence, we conclude that 



log Mb -31o 



n d (b d n d 



71 • • • I 



1 



-H l0 § 
n * * 



13 
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(29) 



+ ^ ¥nd EX n ^R dun j X n -—^j 



-f($)-f a h(»M+o(n-* } y 2 

< 0(n 2 ), 

as n — > oo uniformly over (4>,9i, ■ ■ ■ ,6d) £ Q n .e- Also, we observe from (27) 
that 



(30) <2Var< 



TT d d> d n d 



d \ -i > 

i=l / 



= 0(n~ d ), 

as n — > oo uniformly over /3q,o <4 , ,<f ) < 01,0 and Aj.t < #t < /3i.t , i = 1, • • • , d. 
Finally, since < v < (d-2)/(d + l), it follows from (28), (29) and (30) that 

p{ sup {log[L n (4>, §!,..., 9 d )) ~ log[L„(0, 1} ..., 9 d ))} > log(e)} 

^ 0(n~ d + 2 ) = 0(n( d+1 ^- d + 2 ) -> 
(0,ei,...,e (1 )en n , e 
as n — > oo. This proves Theorem 10. □ 

APPENDIX A 

Lemma A.l. With the notation of Theorem 1, the generating functions 
G T and G T * of {r k : k > 0} and {r^ : k > 0} are, respectively, 

r, , s ^r k z k a + bz ^rf.z k a + bz 



^ u k {l-zf ^ u k {l-zf 



Proof. For all integers k > 0, we have 
J _ = ( w - wu 2 - w 2 u 2 )l{k > 1} + [1 - (1 + w) 2 u 2 ]l{k = 0}, 
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where T{-} denotes the indicator function. Hence 

fc=0 " fc=l fc=l 

which implies that 



1 — (l + -u;) 2 u 2 (w — wu 2 — w 2 u 2 )z 

G T * [Z) = i . 

1 — Z (1 — Z) 1 

Similarly, for all integers k > 0, 

= (w- 2wu 2 + wii 4 )I{fc > 1} + [1 - 2(1 + w)u 2 + (1 + 2w)u 4 ]l{k = 0}. 
Hence 

OO OO fe — 1 00 

fc=0 fc=l fc=l 

which implies that 

(n 2 - l) 2 + 2wu 2 {u 2 - 1) u;(m 2 - l) 2 z 

t(z)_ r^i + (i-z) 2 • □ 

Lemma A. 2. For 1 < m < n — 1, Zei Ag tn - m be the mxm matrix defined 
as in (11). Then with the notation of Theorem 1, 

\A e ,n-, m \ = (-l) m+1 {(««i +u 2 b)[(b + r- 1 u)ai] m - 1 

- (aa 2 + u 2 b)[(b + r„in)a 2 ] m ~ 1 }/(ai - aa)- 

Proof. Using elementary row operations to reduce Ag jn . m to a diagonal 
matrix, we obtain 

m— 1 

|^0,n;m| = (-T-l)" 1 " 1 ^-! + (-T-l) m ~ 2 ^ T m-fci-l^fei-l 

fel=l 

m— 1 fei— 1 

+ (-T_i) m ~ 3 ^-^-lTfe^fca-ir^-l 
fc 1= 2fc 2 =l 

+ ••• 

m— 1 fci— 1 fc m _4 — lfc m _3 — 1 

~ r -l E E " ' E E r m-fei-l7"fci-fc 2 -l "' 
fc 1= m— 2 fc2=m— 3 fc m _3=2 fc m _2=l 

^~km — 3 &m — 2 1 ^~km — 2 1 

+ T m 
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Letting G T {z) be as in Lemma A.l, we observe that |^40,n;m| equals the 
coefficient of z m ~ l in the series expansion of 

(-T.mr^Griz) + (-T_ lU ) m ~ 2 zG 2 T {z) 

+ (-T. lU ) m - 3 z 2 G 3 T {z) + ■■■- T.mz^G™- 1 ^) + z m ~ x G™{z) 

m— 1 

= {-r. lU ) m - l G T {z) (-r^ur k z k G k T (z) 

k=0 

= (-r_ 1 «) G T (*) (_ T _ lti )-i zGr (s) - 1 • 
Following the proof of Theorem 1, we conclude that, for 1 < m < n — 1, 
, _ (-l) m+1 (^ + ^m) m " 1 [Q(ar-a2 1 )+^ 2 &«~ 1 -a™" 1 )] 

«1 — «2 

This proves Lemma A. 2. □ 

Lemma A. 3. Wita the notation of Theorem 4, for 3 <i < j - 2 <n — 2, 
\Re,n;-i,-j;i+l\ 

f_ 2 «f_i[u;(l - 4u 2 + 3u 4 ) + (1 - u 2 ) 2 ] 

T-lT-l ~ T T_2 H 7T— r 

[0 + T^\U)Ol\ — T-\U 

t^-2U 2 [w{2 - 6u 2 + 4n 4 ) + (1 - u 2 ) 2 ] 
(b + r_i«)o!i — r^\u 

MT_ 1 it 2 f- 2 (f_i - MT_ 2 )(1 ~ U 2 ) 2 

[(& + r_i«)o!i — r_iti] 2 

(-l)* _1 (aai +u 2 &) r ., , 
c- i — -[{b + T_iu)ai\ l 6 

Ot\ — Ct2 

f_ 2 «f_i[io(l - 4u 2 + 3n 4 ) + (1 - u 2 ) 2 ] 

T-\T-1 - r r_ 2 H -r— r 

(b + r„iuja2 — r_iu 
f^ 2 u 2 [w(2 - 6u 2 + 4u 4 ) + (1 - u 2 ) 2 ] 

(b + T-\U)OL2 — T^\U 

wr_iu 2 f_2('T-i — iif_2)(l — it 2 ) 2 



+ 



[(6 + r_in)a 2 — r_iu] 



2 



(-l)^ 1 (aa 2 + « 2 6) r ,, , w-3 

^[(6 + r_iu)a 2 ]* • 



ai — «2 
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Proof. First, for 1 = 1,2 and i — 1 < m < n — 1, define 

i-2 

Si sm ,i = ^2[{m-k) - 2(m - fe + l)u 2 + (rn - A; + 2)u 4 ] 



k=2 



x uw^rriCfe + r-iu)"-^ 



k-2k-2 



Then 



Si,™,* - (6 + r_in)a;u r_j S i>ms i 

i-2 



[(m -k)- 2(m - fe + l)u 2 + (m - A; + 2)i 



k=2 



i-2 

- ^ [(m - k) - 2(m -k + l)u 2 + (m - k + 2)u 4 ] 

fc=2 

m— fc— 1 — fc— 1 /» i _ \fc— 1 fc— 1 

x u;?/ r_j [b + T-iuj a, 
[(m - 2) - 2(m - l)u 2 + mu A ]wu m ~ 2 rZf 

- [(m - i + 2) - 2(m - i + 3)u 2 + (m - % + 4)u 4 ] 
x wu m - i+1 TZi +1 (b + T-iuy- 3 ai- 3 

i-2 

- 2(1 - 2u 2 +n 4 )w;n m - fc rr 1 fc (6 + r„ 1 n) fc - 2 af- 2 . 



\fc-2 fc-2 



fc=3 

Hence, for i > 3, 

S'i.m.J _ w[( m -i + 2)- 2(m - i + 3)u 2 + (m - i + 4)u 4 ] 
u m r_iu[(6 + r_in)a/ — r_iu] 

^ + r_in)a/ 1i_3 



+ 



ui[(m — 2) — 2(m — l)u 2 + mu 4 ] 
r_i«[(6 + T-\u)a.i — T-iu] 

w(l-u 2 ) 2 r(6 + r_in)a; ri - 3 



(31) 



[(b + r_in)a; — T-\u} 2 [ T-\U 

w(l — u 2 ) 2 (b + T-iu)ai 
T-iu[(b + T^\u)ai — T-\u] 2 

w[(m -i + 2) - 2(m - i + 3)u 2 + (m-i + 4)ii 4 ] 
T-\u[{b + T-\v)ai — T-iu] 
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2\2 



+ 



w(l — u ) 



[{b + T-\u)ai — T-\uf 



(b + T-\u)ai 



T-\U 



i-3 



w[(m — 2) — 2(m — l)u 2 + rmx 4 ] 
T^\u[{b + T_iu)ai — t_\u] 

w(l — u 2 ) 2 (b + T-\u)ai 
r_iu[(b + T_\u)ai — T_i1i] 2 

Next we observe from Proposition 2 and Case 3 that, for 3 < i < j — 2 < n — 2, 

\Ro,n;-i,-j;i+l\ 

= ( — I) 1 (Ag^ n -_i^j;i + i)i^ + i Y[{Ae : n;-i,-j;i+l)i+l~l,i~l 



1=0 



+ k {Ae in ;-i,-j;i+l)k+l,i+l\Re,n;- 



-i,-j;k\ 



k=l 



(32) 



i-k-l 

x [ (Ae,n;-i,-j;i+l)i+l-l,i-l 
1=0 

= (-1)*(1 + ^)n i+1 f_ 2 f_ 2 rr i 3 

+ (-lf-^l + (i - l)w;]n^ 1 f_ 2 f_ 2 rr i 3 
i-2 

+ y^(-l) t ~ fc Ti-fc|^e,n;-i,-j;fc|''"!-'l fc ~ 2 'r-2^2 
fc=2 

— To\Re,n;-i,-j;i-l\^-2 + ^-1 |^6»,n.;-i,-j;i| 

fc=2 

— TO 1-^61,^-1 1 ^-2 + T-l |-R0,n;-i,-i-l;i| • 

We further observe from Theorem 1 that, for i— 1 <m <n — 1, 

i-2 

,n:k 



k=2 



i-2 



(1 - U 2 ) 2 J2(-ir k U m - k Tl?\R9,n;k\ 
k=2 
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i-2 



+ J2(- i y k i( m ~ k )~ 2 ( m ~ k + 1)« 2 + {m-k + 2)u 4 ] 
k=2 

X WU m - k TZi\Re,n;k\ 



(33) 



u m - l (l-u 2 ) 2 (aai + u 2 b) 
t_i(q!i - a 2 )[ai(b + T-iu) - t~iu] 

u m ~ 1 {l-u 2 ) 2 {aa 2 + u 2 b) 



ai(b + T-iu) 

T-\U 



i-3 



1 



a 2 (b + t-iu) 



+ 



r_i(ai — a2)[a2(b + r_iu) -r_i«] 
(aai + u 2 b)Si t m tl - (aa 2 + u 2 b)S iim>2 



T-\U 



i-3 



ai — a2 

It follows from (32) and (33) that 

\Re,n;-i-j;i+l\ = ("l)*" 1 '^-1^-1 'T- 2 T- 2 - f |i?0, n ;i-l |"?-2 
+ f-l|-R0,n;-i,-(i+l);il 

(-ljVl^f-af-ati^^l - « 2 ) 2 (a«i + n 2 6) 



+ 



(ai - a 2 )[ai(b + T-iu) — t-iu] 
ai(6 + r_iu) V" 3 _ 

(-ljVl^f-af-a^-Hl - n 2 ) 2 (aa 2 + ^ 2 fr) 
(ai - «2) [0:2(6 + r_ 1«) — T_ru] 

a2(6 + r„i^) y" 3 _ 

| (-l)Vri 2 T^2T-2[(Q«i +M 2 &)gj,j,i - (a«2 +^ 2 ^)^,i, 2 ] 

ai — a 2 

and Lemma A. 3 results using (14) after some algebra. □ 

Lemma A. 4. Let be defined as in (17). Then the generating function 
of : k > i + 1} is given by 



fc=i+l 



(1-Z)= 



where 



a 5 



f_ 2 f [w(2 - 6u 2 + 4u 4 ) + (1 - u 2 ) 2 ] 



+ 



r_iu 2 [(6 + r_i«)«i — t_iu] 
jtff_ 2 fo(l - u 2 ) 2 



u[(6 + T-\u)a\ — r_in] : 



+ 



T_lli 4 
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f-gr-iKg - 8u 2 + 5n 4 ) + (1 - tt 2 ) 2 ] 
r_iu[(6 + T-iu)ai — T-iu] 

{aa\ + u 2 b) \{b + t^u)^ 1 ' 



WT-2?-i(1 



,2\2 



[{b + T-\u)ol\ — T-l'uf 



r_!(ai - a 2 ) 



f- 2 T [w(2 - 6u 2 + 4u 4 ) + (1 

T-lU 2 [(b + T_iu)a 2 — T-iu] 



2\2] 



+ 



^T_ 2 f (l 



u 



2\2 



u[(b + T-\u)ot2 — T-iu}- 



+ 



T_i« 4 



f_ 2 f_i[«;(3 -8u 2 + 5u 4 ) + (1 -u 2 ) 2 ] wf_ 2 f_i(l - u 2 ) 2 



T^\u[(b + T_ili)02 — T_in] 

(aa 2 + u 2 6) [(6 + r_iu)a 2 " 11 " 
T^!(ai - a 2 ) I T-iu 



[(b + r_i?i)a 2 — T-iu}' 



wt- 2 u 2 (2 - 3u 2 + u 6 )[w(2 - 6u 2 + 4u 4 ) + (1 - u 2 ) 2 ] 
2r_in 4 [(6 + r_iu)ai — T-\u\ 

w 2 t_ 2 (2 ~ 3u 2 + u 6 )(l - u 2 ) 2 wf (2 - 3u 2 + u 6 ) 

~T~ ~ 77~t ' ~ To 



i/' 



2u[(b + r_in)ai — T-\u\' 

w(l -n 2 ) 2 f_ 2 f_i _ 
r_iti[(6 + r_iii)ai — r_iit] 

(oai + tt 2 6) \(b + r_ 1 n)ai 11 " 



2r_iu 4 
(1 - 3-u 4 + 2u 6 )f_ 



r_!(ai - a 2 ) 



T-\U 



2\2l 



wf_ 2 u 2 (2 - 3u 2 + u 6 )K2 - 6u 2 + 4-u 4 ) + (1 - u 
2r_iu 4 [(6 + T-\u)ct2 — T-iu] 

w 2 f_ 2 (2 - 3u 2 + u 6 )(l - u 2 ) 2 wf (2 - 3u 2 + u 6 ) 
2n[(6 + r_iu)a 2 — r_iti] 2 2r_iu 4 

w(l - u 2 ) 2 f_ 2 f_i w(l - 3u 4 + 2u 6 )f_ 



r_iu[(6 + r_iii)a 2 — r_iu] 
(aa 2 +u 2 6) [(6 + r_iM)a 2 1J 



r_i?r 



r^!(ai - a 2 ) 
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Next, let^k be defined as in (20). Then the generating function Gg of : k > 
i + 1} is given by 

^ h z k Z i+1 (a; + b;z) 



fc=i+l 



where 



{i-zY ' 

f-ifb | f_ 2 f [w(l-4M 2 + 3u 4 ) + (l-u 2 ) 2 ] 
r_iu 4 r_iu 3 [(6 + t_iu)qi — r_i-u] 

wr_if_ 2 f (l - it 2 ) 2 
r_!?i 2 [(6 + r_iu)cti — t_!u] 2 

fif-2 f_ 2 f_ 2 [w(3 - 8n 2 + 5m 4 ) + (1 - m 2 ) 2 ] 
T-iu[(b + T-iv)a\ — r_i«] 

wf_ 2 T_ 2 (l - M 2 ) 2 



[(6 + T-\u)a\ — r_iu] 2 
aai+u 2 6 [(6 + T_iu)ai 1? 



r£ 1 (ai - a 2 ) 



r_i« 



f-ifp | f_ 2 f [w(l-4n 2 + 3M 4 ) + (l-n 2 ) 2 ] 
r_in 4 r_iu 3 [(6 + r_iu)a 2 — r_i«] 

it?r_if_27b(l - u2 ) 2 



+ 



r_iu 2 [(6 + r_iu)a 2 — r_in] 2 

fif_2 r_ 2 f_ 2 [w;(3 - 8u 2 + 5u 4 ) + (1 - u 2 ) 2 } 



T-\u[(b + r_in)a 2 — T-iu] 

I0f_2f_ 2 (l -U 2 ) 2 



a«2 



+ u 2 6 



and 



r£ 1 (o;i - a 2 ) 
iyf_i(2 - 3u 2 + u e 



[(6 + T_iu)a 2 — r_iu] 2 
(6 + r_in)a 2 1 ' 



+ 



+ 



2r_iu 4 

u;f_ 2 (2 - 3u 2 + u 6 )[w(l - Au 2 + 3u 4 ) + (1 - u 2 ) 2 } 
2r_iu 3 [(6 + T-\u)a\ — r_iu] 

w 2 f_ 2 (2 - 3u 2 + u 6 )(l - u 2 ) 2 twf_2f_ 2 (l - u 2 ) 2 



2u 2 [(6 + T-iu)ot\ — T-iu] 2 T-iu[(b + T-iu)ot\ — T-\v\ 
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w{l - 3u A + 2n 6 )f_ 2 

T-lU 3 



au\ + u 2 b 
rL^cti -a 2 ) 



(b + T_iu)a\ 



r_iii 



i-3 



wf_ 1 {2-Zu I + u b ) 
2r_iu 4 

iof_ 2 (2 - 3u 2 + w 6 )Hl - 4« 2 + 3u 4 ) + (1 - u 2 ) 2 ] 



+ 



+ 



2r„iii 3 [(6 + r_iu)a2 — T-itt] 
w; 2 f_ 2 (2 - 3u 2 + u 6 )(l - u 2 ) 2 wr- 2 f- 2 (l ~ u 2 ) 2 



2u 2 [(b + T_ill)a2 — T-lU} 2 T-\u[(b + T-lU)0i2 ~ T-lU] 

w(l -3n 4 + 2n 6 )f_ 2 



T-lU° 



act2 + u b 



(b + T-iu)a,2 



r_iii 



i-3 



Proof. First we observe after some algebra that = £j + i/V +1 and 
a i + b% = ikU~ k — A; > i + 2. Hence, for integers k > i + 1, we 

have 

|| - ^T T { k > i + 2 } = H Z { k = * + 1} + + > t + 2}. 



Hence 



and 



„,k /i.fc— 



fc=i+l 



k=i+2 



— = a^ +l + ( H + b !i ) J2 * 

k=i+2 



Gfc) 



z l+1 (a^ + b^z) 



(1-z) 2 ■ 

Similarly, we observe after some algebra that = £j+i/u' t+1 and a| + b^ 
£,k u ~ k — £,k-i u ~ k+l , k>i + 2. Hence, for integers k > i + 1, we have 

|| - > i + 2} = af{k = i + 1} + (a^ + 6^)X{A > i + 2}. 



Hence 



00 C ~k oo ? Jfc-l 



,k-l = a l z 



t+1 



fe=i+l 



fc=i+2 
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and 



G-Az) 



z i+1 (a e - + bcz 



(1-zf 



This proves Lemma A. 4. □ 



Lemma A. 5. Let < /3 < 0,0 < fii < oo. Then 



n 



^( R ln R eJ 



w 
1 



1 

+ - 

n 



as n — > oo uniformly over (3q<0,0 < (3\. 



2 1 



+ 0(rT 2 ), 



Proof. Let = 0/rt and and m n = \_(n + 3)/2j , where |_ - J denotes the 
greatest integer function. We further define 



e ""r_i 



e u 



1 



and 



ai 



02 



' + r„in)ai «i 2 + -v/3 . 



n 



e"""T_i 
(6 + t„iu)q; 2 



«2 



2- \/3 



l + 0( - 

n 



as n — > oo uniformly over (3q < ,0 < (3\. Using Mathematica, we observe 
from Theorem 6 that, as n — > oo, 

( R e,n R 6,n)n,n = ( R § , n R 8,n) ±A 



i=l 



(34) 



3(70226 + 40545^) 
(664626 + 383722 v^H 



1713 + 989^ 989 + 571^ (0\ 2 



+ 



4(26 + 15^)2 4(26 + 15\/3) 2 \^ 
2702 + 1560\/3 3 



4(26 + 15^)2 \6> 



+ 0(n~ 1 ), 



as n — > oo uniformly over (3q <Q,Q < f3\. Also, we observe from Theorem 7 
that 



( R 6,n R e)i)i,i 



XX 1 + K ~il*) e 

i=i 
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-\i—j\w( d— 1 ' 



+ [1 + (i - 2)«)] e -^- 2 )*(E^) ri _ i+1 , n _ 1 
+ (l + ^) e --(O i _ M + (0 M 



+ (l + ?i ;) e --(i?-i) i)i+1 + [l + ( 



re — ? we 



( R e,n)i,r 



+ 



c 6 V 



Ttxa.{{ai - a 2 ) 



(aa\ + — (aa 2 + ( — 

Oil 



n—i—l 



x{r 4 a (l + 2w)-r 5 a (l + w) 



+ + (< - 3)*] - r l ai + 2)w]}/{<(1 - r Qi n 

C6,2ii 2 r* (aa 2 + « 2 6) 



Tf 1 af(ai — a 2 )(aai + u 2 &) 



(aai + u 2 b) ( — J — (aa 2 + u 2 6) ( — 



n— i— 4-i 



x{r- 1 [l + (i-3^]-r- 2 [l + (i-2)w;] 



-i+3/ 



(l + *)+^ +2 (l+2*)}/{(l-r a2 n 



+ 



C^iu 2 r a l(aai +u 2 b) 

rliaf (ai - a 2 ) 
xK| 2 (l + 2zZ,) 

-^f(l + *)-r™[l + (n-iH 



+ r2+ 1 [l + (n-z-lH}/{(l-r ai n 



g6,m 2 r-;(aQ 2 + n 2 fe) /a 2 x " 
rf 1 af(ai - a 2 ) 



xKt 2 (l + 27i)) 



,i+3 



(l + w) 



rl 2 [l + {n-i)w] 



+ r n a t 1 [l + (n-i-l)w]}/{(l-r a2 ) 2 } 

C Gi2 u 2 r~l(aai + u 2 b)(aa 2 + u 2 b) fa 2 
T"l 1 af(o!i — a 2 )(a«i + u 2 6) 



x{rj+ 2 (l + 2ri3) 



< 3 (1 + *) 
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[1 + (n - i)w] 



+ r n + l [l + {n-i-l)w]}/{{l-r ai T} 



C &)2 u r a l(aa 2 + u 2 b)(aa 2 + u 2 b) fa 2 
T-i<^i {ot\ — ot2){aai + u 2 b) \&i 



n+i— 4 



- rSJl + (n - i)w] + r£ L [l + (n-i- I)w}}/{(1 - r a2 Y} 

as n — > oo uniformly over 5 < i < n — 2 and (3q < 9, 9 < (3\. Consequently, it 
follows from (35) that 

n-2 



H ( R e,n R e)i)i,i 
i=m n 

= (n-m n -l)((^j 



10 

1 



+ 0(n~ 2 )} 



(36) 



3(36810643322 + 21252634831^3) 
(1300170624726 + 750653860178^3)^ 

6322680 + 3650401^3 



128934018 + 74440090^3 

39176289 + 22618441^3 / 9 \ 
128934018 + 74440090^ V Oj \9 



+o(n- 1 ; 



as n — > oo uniformly over (3q<Q,Q<(}\. Furthermore, if n is an odd integer, 
we have 



(37) 



(R 



W 



+ 0(n~ 



as n — > oo uniformly over (3q < 9,0 < Pi. Next we observe from Theorems 
6 and 7 that 



( R 9,n R 6^)n—l,n-l 



= 53(1 + | n _ i _ il^e-l"- 1 -^*^^! 



4G 



(38) 
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3(3351044259 + 1934726305^3) 
(25012534866 + 14440993738^3 jw 

8016837 + 4628523^3 
+ 34547766 + 19946162^3 

^ 3(6584767 + 3801717^) 2 



+ 



34547766 + 19946162^3 \9 
4(1694157 + 978122^3) f9\ 3 



34547766 + 19946162^3 \9 



7T +0(n^), 



as n — > oo uniformly over (3q < 9, 9 < (3\ . Finally, we conclude from (34), (36), 
(37) and (38) that 

1 n 2 n 1 /9\ 3 

-T,( R 8,nRe*nki = ~ E (Re^e'nhi + - ( q ) !{™n = (n + 3)/2} + 0(n" 2 ) 

i=l i=m n 



2(n — m n — 1) 



1 

+ - 

n 

1 



n 



w 
1 



31 « 



+ - U ) I{m„ = (n + 3)/2} + 0(n-'), 
as n — > oo uniformly over (3q <8,9 < (3i. This proves Lemma A. 5. □ 



Lemma A. 6. Lei > and Rg n be as in Section 2. Then as n — > oo, 



and 



1 



tr 



1 

-tr 

n 



3 2(6 + 2) / _. _ 2 . 
-- + ^ - + 0(n 2 ) 



r^ 2 B-Mw 

V d6 2 ' n 



Proof. Since it^ „-Re, n equals the identity matrix, we have via differen- 
tiation 

1 IV <9 „ i \ „ 1 1 



tr 



n 



09 



R e,n ) R e,' 



■tr 



R 



e.n 



(39) 



i=i j=i 



d_ 

89 





Re.n 



'n^X\d9 Re ' n )^ {R ^ k 



h3 
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1 " « (^2^ 
1=1 3=1 



and 



tr 



n 



(40) 



n 



n 



E 

i=l 

E 

i=i 



\d9 2±Le > n 





R, 



-E 

i=i 



Ru 



(96 



#0,n ] Rf) „ 



i.n 



0_ 

00 



Re, 



/cP_ 

\de 2 



Re.n 



where, for all 1 < i,j < n, 



(41) 



9 R 

de Re 



\de 2 



Re, 



1,3 



n 



(i-j) 2 (l-w\i-j\)e-^-i\ w 



n- 



Using (41) and Theorems 6 and 7, Lemma A. 6 is proved by expanding the 
right-hand side of (39) and (40) as a power series in w using Mathematica. 
We refer the reader to [5] for more details. □ 

Lemma A. 7. Let < (3 < 6, 8 < fix < oo. Then 

l^[(R^e,n) 2 ]=0(l), 
as n — > oo uniformly over (3q <0,9</3i. 

Remark. The proof of Lemma A. 7, though conceptually simple, is ex- 
tremely tedious and the symbolic computation software Mathematica fea- 
tures significantly in the evaluation of the error terms. A detailed proof can 
be found in [5]. A much abbreviated proof is given below. 

Proof of Lemma A. 7. First, it is convenient to note that 

( R 6,n R 0,n)h3 = ( R e,n R 9,n) n - i + 1 ,n-j+^^ 

whenever 1 < i, j <n and 

i( R e,n R e,n)\i = [( R e^n R e,n) 2 ^-i+i,n-i+i V 1 < i < n. 
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Also, if n is an odd integer, we have 

[(^n^,n) 2 ](n + l)A(n + l)/2 = (jf + Ofa- 1 ), 

as n —> oo uniformly over (3q < 0,0 < (3±. Hence, writing m n = [(n + 3)/2j , 
we obtain 

^[(i^x-,,) 2 ] 

= ~ E [(Re^n)\i + -{- e ) X{m„ = (n + 3)/2}+0(n- 2 ), 

i=m n 

as n — > oo uniformly over /3q < 0, < Consequently, it follows from 

- E m^,n)%=°o-)> 

i=rn„ 

that 

as n — ► oo uniformly over (5q <0,0 < (3\. □ 

APPENDIX B 

For i = 1, 2, we define 

f_ 2 f [u>(2 - 6u 2 + 4u 4 ) + (1 - u 2 ) 2 } 



Cl,i 



T__ili 2 [(6 + r_iu)aj — T-i-u] 

wf_ 2 r (l -n 2 ) 2 f f - Tjj-i 

u[(b + r_iu)aj — r_iit] 2 r_iu 4 

f_ 2 T_i[«;(3 - 8u 2 + 5u 4 ) + (1 - u 2 ) 2 ] «;f_ 2 f_i(l - u 2 ) 2 



C 



2,i 



r_in[(6 + r_iu)aj — r_iu] [(6 + T_itt)aj — r_iu] 2 ' 

™f_ 2 (2 - 3u 2 + u 6 )[w(2 - 6u 2 + 4u 4 ) + (1 - u 2 ) 2 ] 
2r_iii 2 [(6 + T_iu)aj — T-\u] 

w 2 f- 2 (2 - 3u 2 + u 6 )(l - u 2 ) 2 wf (2 - 3u 2 + u 6 ) 

^ FTt ! \ To 



2n[(6 + T_iu)aj — r_iii] 2 2r_iu 4 

- u 2 ) 2 f_ 2 f_i w(l - 3u 4 + 2u 6 )f_i 



r_iu[(6 + r_iti)aj — T-\u] T-\u 



3 
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f_ 2 u[w{\ -Au 2 + 3u 4 ) + (1 -u 2 ) 2 ] wt^u 2 f^ 2 (l - u 2 ) 2 
C3A = r_i H 71— ; 1- 



C5 ■ 



(6 + T-\u)ai — T-iu [(b + r_i«)ai — r_iu] 2 ' 



c = f_i7b | f- 2 f [w(l - 4n 2 + 3n 4 ) + (1 - n 2 ) 2 ] 
' r_iu 4 r_in 3 [(6 + r_iu)aj — r_i«] 

wr_if_ 2 f (l - n 2 ) 2 
r_iu 2 [(6 + r_iu)aj — r_in] 2 

nf-2 f_ 2 f_ 2 [w(3 - 8u 2 + 5n 4 ) + (1 - u 2 ) 2 } 
r_iu 4 r_!u[(6 + T_iu)ai — r_iu] 

wf_ 2 f_ 2 (l - tt 2 ) 2 
[(6 + r_!n)aj — t^\u] 2 ' 

u;f_i(2-3u 2 + ti 6 ) 
2r_iu 4 

urf_ 2 (2 - 3u 2 + u 6 )[u;(l - 4u 2 + 3u 4 ) + (1 - u 2 ) 2 ] 
2r„iu 3 [(6 + T_iu)aj — r_iu] 

w 2 f_ 2 (2-3u 2 + u 6 )(l-u 2 ) 2 
2u 2 [(fe + r_i«)aj — r_in] 2 

wf_ 2 f_ 2 (l -n 2 ) 2 u;(l - 3u 4 + 2u 6 )f_ 2 



+ 



+ 



T-\u[(b + T_iu)a.i — T_l«] T_iU 3 

f_ 2 uf_i[u;(l - 4u 2 + 3u 4 ) + (1 - u 2 ) 2 ] 



C6,i = f-lf_i — fof_2 + 



+ 



(6 + r_iu)aj — r_iu 

f„ 2 f_ 2 n 2 [w(2 - Gu 2 + 4n 4 ) + (1 - n 2 ) 2 ] 
(b + r_iu)aj — r_iii 

wr_iu 2 f_ 2 (f_i — ur_ 2 )(l — u 2 ) 2 
[(b + r_iu)«j - r_iu] 2 



_ f_ 2 ii 2 [w(2-6u 2 +4ii 4 ) + (l-ii 2 ) 2 ] u;r„if„ 2 u 3 (l -u 2 ) 2 

[(5 + T_in)aj - r_iu] [(6 + r_itt)aj - r_in] 2 °' 
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